Superradiant and Dark Exciton States in an Optical Lattice within a Cavity 
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We study ultracold atoms in a finite size one-dimensional optical lattice prepared in the Mott 
insulator phase and commonly coupled to a single cavity mode. Due to resonance dipole-dipole 
interactions among the atoms, electronic excitations delocalize and form excitons. These exciton 
modes are divided into two groups: antisymmetric modes which decouple from the cavity mode 
forming dark states, and symmetric modes significantly coupled to the cavity mode called bright 
states. In typical setups the lowest and most symmetric exciton is coupled to the cavity photons 
much stronger than the other bright states and dominates the optical properties response of the 
atoms {superradiant state). In the strong coupling regime this superradiant state is coherently 
mixed with the cavity photon to form a doublet of polariton states with the Rabi splitting. 

PACS numbers: 37.10.Jk, 42.50.Pq, 37.30.+i 



The strong coupling regime of a single atom to a sin- 
gle cavity QED mode has been achieved experimentally 
in many different setups [iL and even more extensively 
been studied theoretically As recent breakthrough 
a large Bose-Einstein Condensate (BEC) within a cavity 
has been successfully used to demonstrate strong atom- 
field coupling Experiments now even target to use 
BECs coupled to microwave atom chips [1, Q, where 
strong coupling is expected even for magnetic collective 
coupling on dipole forbidden long lived hyperfine transi- 
tions. Theoretically, in a first approximation the coherent 
coupling of N atoms of a BEC to a single cavity mode, 
with identical single atom-photon coupling /, is well de- 
scribed by the Tavis-Cummings model, which leads to a 
collective coupling of ^/N f [6] . In this model the direct 
electrostatic interactions among the atoms are neglected. 
From the other side, a Bose gas of ultracold atoms loaded 
on an optical lattice realizes the quantum phase transi- 
tion from the superfluid to the Mott insulator Q, which 
is predicted by the Bose-Hubbard model [s*]. Moreover, 
controlling dipole-dipole shifts in an optical lattice clock 
was studied in [§]. These facts bring the fabrication of 
optical lattice ultracold atoms within a cavity very close. 
Such new set-up is of importance for light-matter inter- 
face in quantum information processin g [lOl| . and many- 
particle physics of quantum liquids [Til Il2j|. This system 
exhibits bright and dark collective electronic excitations, 
which promise useful applications for quantum commu- 
nication 13] and quantum memory [l^. 

Optical lattices are formed by prescribed counter prop- 
agating off resonant laser beams interacting with atomic 
motion only in a conservative way. For that atoms in 
such a lattice we consider here only the ground state and 
a single excited state of the atomic electronic excitations. 
We assume that their optical lattice potentials to have 
minima at the same positions and the atoms are local- 
ized in the corresponding first Bloch bands. The optical 
lattice is placed between spherical cavity mirrors, as seen 
in figure (1), so that only a single cavity mode close to 
resonance to the previous atomic transition. The super- 
fluid to Mott insulator quantum phase transition within 



a cavity is studied by us in flS^] , where we confirmed the 
existence of the Mott insulator, but for deeper optical lat- 
tices. Here we concentrate in the case of one-dimensional 
optical lattice in the Mott insulator phase with one atom 
per site, and in the case of low density of excited atoms. 
Different from the Tavis-Cummings model for a BEC 
within a cavity, we include the resonance dipole-dipole 
interaction among the atoms. This interaction gives rise 
to the formation of collective electronic excitations, with 
a superradiant state which is strongly coupled to cav- 
ity photons. Previously we investigated two-dimensional 
optical lattices in planar cavities, and due to resonance 
dipole-dipole interactions in exploiting the lattice sym- 
metry, we got in-plane propagating excitation waves (ex- 
citons) , and within a cavity in the strong coupling regime 
we obtained cavity polaritons [l6j . In the present case for 
finite one-dimensional optical lattice, we have standing 
wave excitons coupled to the fixed cavity mode. 

We present first the cavity photons. Only the lowest 
cavity mode is close to resonance to the atomic transition. 
Hence we consider only the Gaussian beam between the 
mirrors (see figure 1), which is represented by the Hamil- 
tonian He = Ec a^a, where and a are the creation 
and annihilation operators of a cavity photon with en- 
ergy Ec — hi>c, respectively. The electric field operator 
along the waist [l], which is along the optical lattice, is 

E{r, z = Q) = 1^1^ e-^'/< {e a-e* a^), (1) 

where e is the photon polarization unit vector, wo is the 
beam waist, and V is the mode volume, which is given 
hy V = ttwqL/4:, with L the distance between the cavity 
mirrors, as seen in figure (1). Here we neglect the fi- 
nite line width of both the cavity and the excited atomic 
states, to be included later in the linear optical spectra. 

Next we present electronic excitations in the optical 
lattice atoms. An electronic excitation at an atom site 
can transfer among the lattice atoms due to dipole-dipole 
resonance interactions. The excitation Hamiltonian is 

Hex = ^-E'a B^Bn + ^ Jd B^Byn, (2) 
n {n,m) 
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FIG. 1: A cavity of spherical mirrors separated by distance 
L. We consider only the lowest Gaussian beam of waist ujo. 
The optical lattice forms of external lasers is located parallel 
to the waist. 

where and B„ are the creation and annihilation oper- 
ators of an electronic excitation at site n with transition 
energy Ea = hua, respectively. At low number of elec- 
tronic excitations we treat only a single excitation at a 
time. Namely the operators _B„ can be assumed to be- 
have as bosons. We take into account only the interac- 
tion between nearest neighbor sites. The energy transfer 
coupling parameter is Je = ^^^^^^3 (l ^ 3 cos^ (?) , where /i 
is the atomic transition dipole, which makes an angle 
with the lattice direction. For isotropic atoms is fixed 
by the cavity photon polarization direction. 

We consider a finite lattice with N sites. The sites 
are labeled by n = 1, • • • , A^. In order to use fixed 
boundary condition, we add two additional empty sites, 
n = and n = A'' + 1, as in figure (2). Therefore the 
above excitation Hamiltonian can be diagonalize in us- 
ing the transformation i3„ = y^^^X^fe ^i'^ (iTTT^) 
where the collective excitation modes are labeled by 
k = 1, - ■ ■ ,N, which are considered as standing excitons, 
and are presented in figure (2). The diagonal Hamilto- 
nian reads H^^. = '^^.Ef^ Bj^Bk, with the energy disper- 
sion Ek = Ea + 2Je cos (^-^^^ ■ In place of discrete atom 
levels we get an energy band of band width 4J0. 



k=3 k=2 k=l 




FIG. 2: One-dimensional optical lattice of one atom per site 
Mott insulator, with lattice constant a. In using fixed bound- 
ary condition we get standing excitons. The first three exciton 
modes axe plotted. 

For the coupling between the electronic excitations 
and the cavity photons we use the dipole interaction 



Hi = —jl-E. The dipole operator is defined by /t = 
P-Tlm i^n + ^n)j ^Lud wc apply the rotating wave ap- 
proximation. As the photon polarization at the waist 
is in the optical lattice direction, we choose here 9 = 
and then n = {e- p). The maximum of the Gaussian 
is taken in the middle of the optical lattice. Then the 
atom positions r„ = na — L/2. with the lattice length 
L = {N + l)a. For wq enough larger than L, we use 
e,~^n/'^o ~ 1. We need also the summation over n, 
which is I]« sin ^-j^fc^ = cot (^ 2{n+i) ) ^"'^ 
(fc = 1, 3, • • • ), and zeros for even fc-s, (fc = 2, 4, • • • ). Fi- 
nally we get Hi = J2k{odd) {fk + f 'l a)B,^ , where 

= ~^\J eo^iV-fTJ "^ot (2(JV+1))- 

The excitonic states are divided into two parts. The 
antisymmetric states of even fc-s, (fc = 2, 4, • • • ), with odd 
number of nodes. These states are decouple to the cavity 
photons, and are considered as dark states. The symmet- 
ric states of odd fc-s, (fc = 1, 3, • • • ), with even number of 
nodes. These states arc coupled to cavity photons, and 
are considered as bright states. Here we show that the 
first mode (fc = 1), which has no nodes, is a superradiant 
state, and the rest states (fc = 3, 5, • • • ), are weakly cou- 
pled to the light. The things are clear due to the function 

cot ^ 2(jv+i) ) which decays very fast for small fc-s, where 
|/i|/|/fe/i| = fc- The oscillator strength is proportional 
to l/fep, hence the oscillator strength of the first mode 
(fc = 1) is stronger by fc^ from the (fc ^ 1) state, for small 
fc-s. For example the oscillator strength of the first state 
is stronger 9 times from the third state. As the oscilla- 
tor strengths are summed to one, we conclude that, for 
{N ^ 1), the first state (fc = 1) includes 0.81 of the sum, 
and the second bright state (fc = 3) includes 0.09. 

To present the above results we consider an optical 
lattice within a cavity for the following numbers. The 

lattice constant is a = 10~^ V^]: the beam waist is wg = 
3 X 10~"^ [M], the distance between the mirrors is L = 
1.5 X 10~^ [M], and the cavity mode volume is F = 
lO""'^" [M^]. The atomic transition dipole is /i = 5 x 
-|^Q-29 [(7_;vf]^ the angle between the dipole and the optical 
lattice direction at the waist is taken to be ^ = 0, and 
which results in the energy transfer parameter of Jq /h = 
—6.8 X 10^ \Hz]. The atomic transition frequency is Va 
■ X 10^^ [Hz]. For the case of = 10^ sites, in figure 
(3a) we plot the shifted exciton dispersion Eh — Ea as 
a function of fc. In figure (3b) we plot the square of the 
exciton-photon coupling |/fep as a function of fc. Here the 
cavity photon frequency is taken to be at resonance with 
the first exciton mode, that is Ec = E\. It is clear how 
the coupling decays very fast for large fc. For the coupling 
of the first exciton we have \fi\/h = 2.55 x 10*" [Hz], and 
for the third exciton we have ] fs\/h = 8.5 x 10^ [Hz]. 

In the light of the above discussion, we can neglect 
the coupling of the cavity photons to the exciton states 
except the coupling to the superradiant state. Therefore, 
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FIG. 3: (a) The shifted exciton dispersion {Ek — Ea)/h vs. fc, 
for A'^ = 10^ and 6 — Q. The dispersion is a cos function with 
band width 4Jo. (b) The squared exciton- photon coupling 
\fk\'^ /h vs. k, for — . It is seen that the first mode has 
the dominant coupling, which is nine times larger than the 
third mode. The even modes are dark. 



the coupled photon-exciton Hamiltonian is 

H = Eex B^B + Ec.a^a + f B^a + f* a^B, (3) 

where we dropped the exciton index. The coupling pa- 
rameter is 



eaV{N + l) 



cot 



2{N + l] 



(4) 



El A\Ar, with 



+ £;e.)/2±A, 
with the detuning (5 = (E^ — 



and the exciton energy is E^^ = Ea + 2Jg cos (^j^^ 

In the strong coupling regime, where the exciton-photon 
coupling is larger than their line widths, the Hamiltonian 
can be easily diagonalized in using the transformation 
A± = X± B -h y± a, to get H = J2, 
the two polariton eigenenergies E^ = {E, 
where A = V'^Tj/p 
Eex)/'^- The eigenstates are coherent superpositions of 
the exciton and the cavity photon. The exciton ampli- 
tudes are X"^ — ±y/(A =f (5)/2A, and the cavity photon 
amplitudes are = //^2A(A =F S). In figure (4a) we 
plot the shifted polariton dispersion E^ — E^x as a func- 
tion of the detuning 6, in using the previous numbers. In 
figures (4b) we plot the excitonic and photonic weights 
for the lower and upper branches. 

To observe the system we need to couple the cavity 
mode to the external radiation field. We do this by con- 
sidering non-perfect cavity mirrors. We introduce cavity 
mode damping rates for the upper and lower mirrors by 
7. As in experiments the cavity is open, we also include 
the cavity photon damping rate Fc directly into the free 
space. The excited atoms have damping rate of Fq, which 
is included here phenomenologically. The external field 
serves us with two thin gs: damping of the cavity mode, 
and input-output fields lu\. 

We interest in the linear optical spectra. For an inci- 
dent given field with a fixed polarization from the upper 
side of the cavity and normal to the mirror, in figure 
(5a) we plot the transmission spectra, and in figure (5b) 
the refiection spectra. We used the damping frequencies 
Fq = Fc = 7 = 10^ [Hz]. The two peaks and the two 
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FIG. 4: (a) The shifted polariton dispersions {E^ — Eex) /h vs. 
the detuning S/h. The vacuum Rabi splitting is obtained at 
5 — 0. (b) The excitonic and photonic weights, I^^P, 
vs. the detuning S/h, for the lower and upper branches. At 
S — the branches are half exciton-half photon. For large 
positive detuning the lower branch becomes exciton and the 
upper becomes photon, and vice versa for negative detuning. 



dips correspond to the two polariton states. We assumed 
here zero detuning, that is (5 = 0, and 9 — 0. 
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FIG. 5; (a) The transmission spectra, (b) The reflection spec- 
tra. The two peaks and two dips correspond to the polariton 
doublet, at 5 = 0. The spectrum line widths result of the 
cavity photon and excited atom line widths. 



For comparison needs, we consider the same previous 
system, but now in neglecting the dipole-dipole inter- 
actions among the atoms. As we neglect also the spa- 
tial change of the cavity mode along the optical lattice, 
this makes the system similar to trapped ultraclod boson 
atoms within the previous cavity, in the middle between 
the mirrors, where the atoms are taken to be independent 
by neglecting the dipole-dipole interactions 3. 

The electronic excitation is described now by the 
Hamiltonian H^x = X)n B^Bn, as before, for low 
excitations the operators -B„ are taken to be bosonic. 
The excitation-photon coupling Hamiltonian reads Hj = 
(/ ^n'^ + /* '^^Bn) , whcrc the coupling parameter is 

xhe total Hamiltonian is given by 



/ 



2eoV 



H = Y,Ea BlBn + E, a^a + Y.{f B^a + f* a) B„ 

n n 

We define the collective excitation operator Bn 
^B,B = ^J:^B„, to get 



(5) 



H ^ Ea B^B + Ec a^a + f B^a + f* a^B, (6) 
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The Hamiltonian 



is similar 



where f = —1,-1/ — tt- 
to the one for the coupled superradiant states and pho- 
tons with the differences in the energies, Ea and 
Eex^ and the coupling parameters. The energy shift 

is Ea — Eex — — 2Je cos (^]vTt) ^'^'^ equal to about 
f.35 X 10*^ [Hz] for N = 10^ and 6* = 0, which is a 
significant shift. The diagonalization of the Hamiltonian 
is exactly as before, and also the linear spectra calcu- 
lations. For the case of resonance Ec = Ea, we have 
\f\/h = 2.8 X 10^ \Hz\. In figure (6) we compare between 
the vacuum Rabi splitting frequency, Oq — 2\f\/h, as a 
function of the atom number for the two cases. It is seen 
that for large atom number the vacuum Rabi splitting 
for independent atoms is larger than that of interacting 
atoms. The difference is about 5 x 10^ [Hz] for N = 10^ 
atoms. We conclude that the vacuum Rabi splitting is 
reduced by the dipole-dipole interactions. 
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FIG. 6: The vacuum Rabi splitting frequency Q.o vs. the atom 
number N, for interacting (full line), and non-interacting 
(dashed line), optical lattice ultracold atoms. The difference 
increase with the atom number. 



and which is N and dependent. In figure (7a) we plot 
the generalized Rabi splitting as a function of for the 
two cases, and with TV = 10"^. It is clear that around 
9 — 54.74° the generalized Rabi splitting for interacting 
case is lower than the non-interacting case, where around 
9 = 54.74° the dipole-dipole interaction is zero, that is 
Jg K, 0. In figure (7b) we plot the generalized Rabi split- 
ting as a function of N for the two cases, where for the in- 
teracting case we plot for the angles 6 — 0°, 54.74°, 90°. 

In summary we explicitly calculated the coupled atom- 
field eigenmodes of a quantized cavity mode strongly cou- 
pled to collective excitations (excitons) of cold atoms in 
an optical lattice. Including resonance dipole-dipole in- 
teractions strongly influence the system eigenmodes and 
the electronic excitations form collective dark and bright 
states. Within a typical cavity only a single superradi- 
.xio" 
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FIG. 7: (a) The generalized Rabi splitting frequency SI vs. 6, 
for interacting case (full line), and the non-interacting case 
(dashed line), with = 10^. The non-interacting case is 6 
independent, (b) The generalized Rabi splitting frequency 
Q. vs. A^, for interacting case (full lines) at the angles 6 = 
0°, 54.74°, 90°, and the non- interacting case (dashed line). 
Around 6 — 54.74° the interacting case splitting is smaller 
than the non-interacting one. 



For more comparison we plot the generalized Rabi 
splitting for the two cases. We consider the case with 
Ea ~ Ec- Then for atoms without dipole-dipole inter- 
action we have S — 0, and then = 2[f[/h, which is 
a square root function of N, and 9 independent. For 
the case with dipole-dipole interactions, we consider only 
the superradiant state, to get the generalized Rabi split- 
ting, which is = 2A/h, where A defined previously. 



ant mode dominates the linear optical response forming a 
doublet of eigenstates. Note that analogous results also 
hold for any chain of electromagnetic active materials, 
e.g. a lattice of quantum dots within a cavity. 
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